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For compact U(l) lattice gauge theory (LGT) we have performed a finite size scaling analysis 
on N T N^ lattices for N T fixed and N s — > °°, approaching the phase transition from the confined 
phase. For N T = 4, 5 and 6 our data contradict the expected scenario that this transition is either 
first order or in the universality class of the 3d XY model. If there are no conceptional flaws in 
applying the argument that the Gaussian fixed point in 3d is unstable to our systems, estimates of 
the critical exponents a/v, 7/ V, (1 — j3 ) / V and 2 — 77 indicate the existence of a new, non-trivial 
renormalization group fixed point for second order phase transitions in 3d. Such a fixed point 
would be of importance for renormalization group theory and statistical physics. 
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1. Model and Simulation Statistics 

This paper summarizes and extends work of Ref. §\§ about pure U(l) LGT with the Wilson 
action. Multicanonical simulations [^] on 

NiN* lattices, N = N r N^ (1.1) 

are performed with N t = fixed, N s — ► °° or N t = N s . The multicanonical parameters were de- 
termined using a modification of the Wang-Landau (WL) recursion [ph. A speed up by a factor 
of about three was achieved by implementing the biased Metropolis-Heatbath algorithm [Q] for 
the updating instead of relying on the usual Metropolis procedure. Additional overrelaxation [j|] 
sweeps are presently used for simulations on larger lattices. 

Our statistics is summarized in table [T]. The lattice sizes are collected in the first and second 
column. The third column contains the number of sweeps spent on the WL recursion for the multi- 
canonical parameters. Typically the parameters are frozen after reaching e = e 1 / 20 . Exceptions are 
the simulations for which the values are marked by * for which e = e 1 / 22 was used (technical de- 
tails will be published elsewhere). Column four lists our production statistics from simulations with 
fixed multicanonical weights. Error bars as shown in figures will be calculated using jackknife bins 
(e.g., chapter 2.7 of [Q]) with their number given by the first value in column four, while the second 
value was also used for the number of equilibrium sweeps (without measurements) performed after 
the recursion. Columns five and six give the j3 values between which our Markov process cycled. 
Adapting the definition of chapter 5.1 of [§] one cycle takes the process from the configuration 
space region at j8 m j n to j8 max and back. Each run was repeated once more, where after the first run 
the multicanonical parameters were estimated from its statistics. Columns seven and eight give the 
number of cycling events recorded during these runs. We are still producing on lattices with larger 
N s values than those listed in table [j] and are creating entire new series of N s values for N r = 2, 3, 
8 and 10. We perform finite size scaling (FSS) calculations for the critical exponents of U(l) LGT. 
For a review of FSS methods and scaling relations see The observables, which we calculated, 
and their FSS behavior are introduced in the following. 

For the Specific heat one has 

C(/3) = ^[(S 2 )-(S) 2 ] ma ^N^ v . (1.2) 

Polyakov loops Pj are the products of U(l) gauge matrices along straight lines in the N t direction. 
The lattice average of Polyakov loops is defined by 

P = Y f p x - (1-3) 

X 

The maxima of the susceptibility of the absolute value \P\ (called Polyakov loop susceptibility 
henceforth) scale like 

Xmax = ^[(l^| 2 )-(l^l) 2 ] max ~^ /V , (1-4) 

while for the maxima of the derivatives the scaling behavior 

1 d 



a = J_^_/|p| 



N^ /V . (1.5) 
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Table 1: Statistics. 
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Figure 1: Maxima of the specific heat, 
holds. Maxima of structure factors (see, e.g., Ref. []8]]) 
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Y.P(r) exp(ikr) 

r I 

F(k) max -ivf"" =N? /V . 



a = 1 lattice constant) 



(1.6) 



(1.7) 



2. Estimates of Critical Exponents 

The lattice size dependence of the maxima of the specific heat is shown in Fig. [jj For our 
symmetric lattices a fit to the first-order transition form 

C max (N s )/(6N) = co + ay/N + ^/N 2 (2.1) 

yields the estimate with Co ~ 0.00020. Statistical error bars of our estimates can be found in |jl|]. 
Our Co value is 10% higher than the one reported in Rev. [Bp, where lattices up to size 18 3 were 
used for this fit. 



and (1.7) 



In contrast to symmetric lattices we find for N x = 4, 5 and 6 from Eqs. (1.2), (1.4), (1.5) 

l-P 



v 



^2, 



1.5 and 2-tj«2. 



(2.2) 
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Figure 2: Polyakov loop susceptibilities for N z = 4. 



The hyperscaling relation 2 — a = dv with d = 3 implies then aava 0.5. For the Polyakov 
loop susceptibility and N t = 4 Fig. ||| shows our canonically reweighted data, implying y/v = 2 
and, hence, y rs 1. For the reweighting the logarithmic coding from chapter 5.1.5 of ^ was used. 
From the derivative of (\P\) we get similarly j8 « 0.25, which is consistent with the scaling relation 
a + 2/3 + y = 2. While from the structure factor we find 2 — rj consistently with the scaling relation 
y/v = 2-r\. 

The exponents listed are the 3d Gaussian values. In view of statistical errors and expected 
systematic errors due to our limited lattice sizes, close-by values are as well possible. The problem 
with the Gaussian values is that the Gaussian renormalization group fixed point in 3d has two 
relevant operators [10]. So it is unstable and one does not understand why the effective spin system 
should care to converge into it [11]. Therefore, the interesting scenario of a new non-trivial fixed 
point with exponents accidentally close to 3d Gaussian arises. An illustration, which is consistent 
with the data, are the values v = 0.482, a = 0.554, y = 0.94, j8 = 0.253, T] = 0.05. Of course, 
we cannot entirely exclude that finite lattices mislead us and that really large systems turn either 
towards a first order transition or the 3d XY fixed point. 

U(l) LGT in the N t N^,N s > N t geometry has previously only been studied by Vettorazzo and 
de Forcrand [Oj. While their result for N t = 4 are consistent with a second order transition, they 
claim first order transitions for N x = 6 and 8. However, their N x = 6 and 8 data on very large 
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32 




Figure 3: Rescaled maxima of Polyakov loops susceptibilities. 

lattices may simply exhibit the critical slowing down, which is typical for second order transitions 
[jl|]. For N t = 4, 5 and 6 our data support, independently of N r , the same critical exponent. We 
illustrate this here by rescaling the maxima of our Polyakov loop susceptibilities with a common 
factor, so that they become equal to 1 on symmetric lattices. On a log-log scale the results are 
then plotted in Fig. [3] against N s /N t . The behavior is certainly consistent with assuming a common 
critical exponent for all of them (parallel lines are expected for large N s /N t ). The figure includes 
also some preliminary N t = S data, which so far blend in. 



3. Summary and Conclusions 

The litmus test for identifying a second order phase transitions is that one is able to calculate its 
critical exponents unambiguously. With MC calculations and FSS methods this finds its limitations 
through the lattice sizes, which fit into the computer and can be simulated in a reasonable time. 
Our lattice sizes are not small on the scale of typical numerical work on U(l) LGT, for instance the 
lattices used for the cq estimate of [g]. However, in view of the fact that our data do not support the 
generally expected scenarios, it would be desirable to extent the present analysis to lattices up to, 
say, size 32 4 . Using supercomputers this appears feasible. With mass spectrum methods [13], see 
JT4[ ] for recent work, one may investigate the scaling behavior of the model from a different angle. 
Finally, renormalization group theory could contribute to clarifying the issues raised by our data. 
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